When Do Superfluidity and Long Range Order Imply Entanglement? 
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We investigate tacitly assumed relationships between the concepts of super-fluidity (- 
conductivity), long range order and entanglement. We prove that the three are by no means 
equivalent, but that notwithstanding, some rigorous implication can be established between them. 
This leads to three different, albeit frequently related, notions of "criticality" , all of which are ex- 
emplified within the Hubbard model in the low density regime. We use Peierls' method of twisted 
Hamiltonians to link the existence of entanglement to superfluidity and (quasi)-long range order. 
As an application of our formalism, we show that recent experiments with cold atoms already prove 
the existence of the field theoretic, spatial entanglement in two dimensions. More interestingly, 
the appearance of entanglement in these experiments seems to be intimately related to the phase 
transition of the Kosterlitz Thouless type. 

PACS numbers: 03.75.Gg, 03.67.Mn 



INTRODUCTION 



Phase transitions play a very important role in physics because they challenge us to describe these inherently 
macroscopic phenomena in terms of the underlying, microscopic, mechanisms \\. We have many different types of 
! phase transitions, ranging from a solid to liquid to gas transition of large numbers of atoms, to the quantum mechanical 
■ electron (Cooper) pairing signifying the onset of superconductivity. All different phase transitions, however, have one 
' crucial aspect in common. Simply put, by varying external conditions a certain same lump of matter (or light) can 
fundamentally change its macroscopic behaviour. Ehrenfest was the first to classify phase transitions according to 
which macroscopic property in particular suffers this abrupt change. If the first derivatives of the free energy (such as 
the internal energy, or magnetisation) are the observables in question, then the phase transition is of the first order; if 
the second derivatives become discontinuous (e.g. magnetic susceptibility or heat capacity), then the phase transition 
is of the second order, and so on. In spite of this advance, a coherent account of phase transitions was still amiss until 
y \ the application of quantum mechanics to many body physics. 

The first theory of superfluid critical phenomena was developed by Landau [2| and also goes under the name of 
mean-field theory. It suggests that any phase transition can be viewed as an order to disorder transition. To quantify 
it we therefore need to identify a function (defined at all points in space), called an order parameter, whose non- 
zero value identifies that we are below the critical point (e.g. temperature) where order exists. Above the critical 
temperature the order parameter disappears. The existence of a function that stretches across the whole physical 
system in the ordered phase immediately suggests that different parts of the system, even if far apart, are in some 
sense correlated through this order parameter function. This gives rise to the notion of long range order (LRO), 
introduced by Penrose and Onsegar [3| , in order to understand Bose condensation. The condition for the existence of 



LRO is given by 



lim (*t(x)*(x')) => const > (1) 



and for Penrose and Onsegar is identical to the system being in a Bose condensate. The constant to which the two 
point correlations tend in the large size of the system limit is, in fact, the order parameter discussed before. For us 
here the interesting question will be if this long range order implies the existance of entanglement and if so, what 
kind of entanglement might that be. Before that let us first finalise our brief exposition of phase transitions. 

By the time low temperature superconductivity was properly described by Bardeen, Cooper and Schriefer it 
was clear that the LRO does not quite apply to it immediately. Electrons in a superconductor are not correlated 
across such large distances that they cover the whole of the superconducting sample. Yang was the first to realise that 
we need a new notion of order which he termed the "off-diagonal long range order" (ODLRO) [5]. If instead of ^(x) 
representing a single electron at position x, this operator represented a Cooper pair at the same x, ^>(x) = ij;^(x)ipi(x) 
then the Penrose Onsegar condition for LRO would still be satisfied by this new operator. The proper order parameter 
is therefore the wave function for a Cooper pair - a pair of electrons in a spin singlet state. We can then think of 
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a superconductor as a condensate of Cooper pairs, though electrons are still fermions and behave as such. (In BCS 
theory of superconductivity, we can likewise think of the energy gap between the ground state with Cooper pairs and 
the excited states where they start to break up as a good order parameter. The gap disappears exactly above the 
critical temperature, when the superconductor transforms into an ordinary conductor. Both of these formulations are 
equivalent and lead to the same calcualtion of the critical temperature for superconductivity.) 

The next big surprise in understanding phase transition came in the seventies when it was discovered that we 
can actually have a phase transition without the existence of any long range order. This discover was, interestingly, 
predated by a number of different proofs of the fact that continuous phase transitions cannot exist in 1 and 2 spatial 
dimensions @, 0, @| • The reason for this is very simple and it lies in the fact that any long range order is destroyed 
by thermal fluctuations in low dimensions. However, what can happen is that a kind of short range order emerges 
in two dimensions (we will refer to this a quasi long range order) out of a disordered phase. This means that below 
some temperature the correlations exhibit a polynomial drop with the distance like so limi x _ x 'i_, o0 (^ , ^(a;)^ r (a;')) — > 
l/\x — x'\ p , where p is some power. Above this temperature the drop is typically exponential. The mechanism for this 
transition was explained by Berezinskii and Kosterlitz and Thouless [9] (BKT) and has been experimentally verified 
a number of times since. In summary, therefore, we can have long range order, off diagonal long range order as well 
as quasi long range order present in matter and all linked to phase transitions. 

The whole story about order-disorder transitions can and has been applied to quantum phase transitions as well 
[lfl | . These transitions occur at zero temperature (in practice, at low temperatures) and are driven by changes of some 
parameter other than the temperature, such as the external magnetic field in a spin chain or a the doping parameter 
in a high temperature superconductor. Quantum phases can also exhibit long range order and quasi long range order. 
The point here, of course, is that any correlations are now likely to imply some type of entanglement, since we have an 
overall pure state In much the same way that different correlations in three dimensions lead to different crystals, 
we might expect that different types of entangled state will lead to different phenomena in quantum critical regions. 



INDEPENDENCE OF ENTANGLEMENT AND QUANTUM ORDER 

In this paper we address the issue of if and how quantum criticality and entanglement are related. We will find 
that the relationship is not entirely straightforward, but that some concrete conclusions can nevertheless be drawn. 
Before we start our more general discussion, we note that entanglement and quantum ODL are seemingly not strictly 
related to one another: each can exist without the presence of the other. For example, a product of coherent states 
of equal amplitude but at different positions x, Y[ x \ a ( x ))i is spatially a disentangled state (by construction), but it 
exhibits ODLRO (equal to \a\ 2 ). On the other hand, the state of the GHZ type 1 000. .0) + 1 1 1 1 . . 1) is clearly entangled, 
but LRO does not exist. The latter can be seen by the fact that the operator combination ip'ip vanishes for states 
with the same number of excitations such as 1 00) and |11) (Here it is more appropriate to introduce Pauli raising and 
lowering operators, but we leave this for later). Therefore, any kind of criticality that is indicated by LRO can be 
very different from the criticality signified by the vanishing of entanglement. Similar statements can be made about 
entanglement and quasi long range order. Entanglement is, in fact, just a more complex form of quantum order, and, 
just like correlations, it can exist between any collection of spins, be they close to each other or on different sides of 
the investigated system (for a recent review of entanglement in many-body systems see 

If, as above, we can show that entanglement and LRO are independent, why has there been so much work on 
entanglement and phase transitions? The simple answer is that states do not encode all the information about the 
actual physics of the system. The Hamiltonian plays an equally important role. It can easily happen that the states 
we discussed above (coherent and GHZ) are never the eigenstates of the relevant Hamiltonian. And this is exactly 
what will happen in the situations analysed below. We will show that entanglement is, in spite of the above arguments 
to the contrary, intimately connected to quantum order and critical phenomena. 

Our approach here will entirely be based on physically observable effects. The notions of different order parameters 
are useful precisely because they are ultimately linked to some underlying physics. Long range order in supercon- 
ductors, for example, has been shown to imply the Meissner effect as well as the notion of flux quantisation [l3| (in 
superfluids it leads to irrotational flow and quantized vortices [HI)- It is in this sense that it would be desirable to 
view entanglement: can we say that superconductivity (-fluidity), or some other critical phenomena, are in any way 
dependent on entanglement? Here we analyse this question in detail and show that the answer is affirmative. Our 
example system will be very simple, but it contains all the necessary elements to draw some general conclusions. We 
will analyse both classical (i.e. T > 0) and quantum (T = 0) criticality and show that recent experiments with cold 
atoms already confirm the existence of field theoretic, particle number, entanglement. 
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ANALYSIS OF SPINLESS HUBBARD, XX MODEL 

In order to remain close to physics, we consider the so called Bose -Hubbard model. This model is used to 
investigate various superconducting and superfluid behaviours, ranging from High T c superconductors to cold atom 
gases in optical lattices. The model has quite a complex phase diagram [l5j |. but we will only be interested in some 
special regimes relevant to our study of entanglement. The Bose-Hubbard Hamiltonian is given by 

H BH = J^2( b i b ^ + b 4+i) + \ E( n « - X K > (2) 

i i 

where the first term describes the nearest- neighbour site hopping of bosons, (b^ and bi+i are the usual bosonic raising 
and lowering operators respectively), and the second term is the on site repulsion between bosons. In the limit of low 
density, the nf term can be ignored and we obtain the so called XX Hamiltonian: 

H = -j£ « +1 + « +1 - M ]T a\ (3) 

i i 

where a x = b + tf, a v = i(6— 6 f ) and a z = 1-26*6 are the usual Pauli matrices and we can think of \x as the chemical 
potential (iq| . They represent two level systems where the states |0) and |1) are the boson occupation numbers at 
each site (an empty site and one boson in the site respectively [16]). Physically this should be clear, since by making 
the density very low, we preclude more than one particle from occupying each of the sites. Note that this Hamiltonian 
would also be obtained if we investigate the spinless fermion Hubbard model since fermions obey the Pauli exclusion 
principle and therefore we cannot have more than one electron per site. 

The XX model has been extensively studied and its spectrum is well understood through the Jordan- Wigner 
transformation [10]. However, since we would like to understand superfluidity in the original Hubbard model we need 
to know the response of the system with such a Hamiltonian to introducing external perturbations (here we study a 
one dimensional model; two dimensions are needed for proper superfulidity, but the one dimensional discussion leads 
to the same conclusion - see the Appendix for a discussion of this point). In order to study criticality in the XX 
model we introduce the so called "twisted" Hamiltonian [l7. 18|: 



He = -J^eV^ +<r i -<r+ 1 e- tf (4) 



e 



obtained by imposing Peierls phases on each Pauli raising and lowering operator: a — > e l0 a and a + 
We can think of phases arising from an imposition of an external field (like, for example, the Aharonov-Bohm effect, 

where the charge encirculating the vector potential, A, gains the phase e~ l § Adl ). The twisted Hamiltonian therefore 
represents the response of the system to an external disturbance and this is what gives rise to the physics of various 
critical phenomena. The superfluid density can now be obtained from the assumption that the difference between the 
normal and "twisted" energy is - to the lowest order - given by the superfluid kinetic energy. So, when we externally 
weakly perturb the system, its response is to move its superfluid component (if there is such a component). If the 
response is null to the lowest order, our material is not in its critical phase. Translated into mathematics our statement 
reads [lj|: 

(Hg) - (H) = ^Nf s mv 2 (5) 

where N is the number of particles in the superfluid, m is the mass of each particle and v s the superfluid velocity. The 
superfluid density is denoted by f s . The superfluid velocity is given by v s — V0, where 9 is the phase of the superfluid 
macroscopic wave function - the order parameter (and is the same angle that enters the twisted Hamiltonian). We 
will assume that the superfluid phase varies linearly with distance, 9(x) = 9tx/ L. From this definition it is clear that 
the superfluid fraction, f s , can be expressed as: 

_ 2m L 2 (Hg) - (H) 1 (Hg) - (H) 
Js h 2 N 9 2 T JN 9 2 1 1 

which is essentially the energy difference between the perturbed and the original Hamiltonian, divided by the phase 
squared. Note that here we used the fact that J = h 2 /2ma 2 and 9 = 9t/N, where a is the typical lattice spacing. 
Computing (Hg) to the second order in 9, we obtain the superfluid fraction XX to be: 
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The second term, usually known as the superfluid current [18j, vanishes for our translationally invariant Hamiltonian 
(as can easily be checked) . The subscript "2" indicates that the average is a more complicated second order one, but 
the exact details do not concern us here (see (l8j|). The condition for non-zero superfluid density, f a > 0, now leads to 

^(E«+i+^ + i)^° ( g ) 

i 

So the supercurrent and therefore criticality exist as long as the above expectation remains finite, which in turn 
implies that the condition for criticality is [i < J (so long as we are at T = 0). When, on the other hand, fx > J, 
the superfluid density identically vanishes. The critical point is therefore fj, = J and the state above critical point is 
|00...0), signifying that all sites in he system are empty. Physically this is because the chemical potential, /i, is too 
high for any particles to exist (note that the \i influences the form of the ground state and, therefore, it also affects 
the expectation value determining the superfluid fraction, even tough it does not appear explicitly in the equation eq. 

©)• 

Let us turn to the existence of entanglement. Since at zero temperature the state is pure it is easy to determine 
when it is entangled and when not. And exactly the same critical point (as for superfluidity) is identified as the 
borderline, J = fi. We can quantify entanglement in a number of ways, and single site entropy is just one of them. 
This is given by S = — e+ lne + — e_ Lne_, where e± = (1 ± (<j z ))/2. For \i > J, (a z ) = 1, and S vanishes. Otherwise, 
entanglement always exists and peaks at the point when fi = 0. At zero temperature, therefore, superfluidity and 
particle number entanglement [20] are fully equivalent (this now formally proves partially supported claims in [2l|). 
Long range order, on the other hand, vanishes, but short range order does exist as discussed above. The just described 
phase transition at J = [i and T — is an example of a quantum BKT phase transition [Io| ; above criticality the 
long range order vanishes, whereas below it suffers a polynomial drop. 

At finite temperature we can no longer assume that entanglement is quantified by S since the overall state is mixed 
[22I ] . We thus have to resort to the technique of entanglement witnessing. We have addressed this problem in several 



previous publications [23|, [24j, [25| (see also [26|). The result, when applied to the current scenario says that f s > 1/2 
implies the existence of entanglement (this is because the correlations leading to this value of the superfluid density 
can only possibly be explained by assuming entangled states). Below this value of f s we still have the superfluid, but 
it is not clear what we can conclude about entanglement. The two phenomena are therefore at finite temperature 
most likely not fully equivalent, though they still are very much related. At high temperature the condition for 



entanglement can roughly be expressed as /i 2 + T 2 < J 2 |24(. Therefore, at temperatures where superfluidity exists, 
we can conclude that entanglement also exists. 

Nearest neighbour entanglement can likewise be calculated, since all this requires are the knowledge of (of of, -l), 
(of of +1 ) and (erf). The concurrence, C, between any nearest neighbours is then given by C — max{0, |(of of +1 )| — 

J(l + (of of +1 )) 2 - 4(crf ) 2 }. Furthermore, (of of +1 ) = {U-(J,M)/2J Q so that at fj. = T = we have that C = 0.04 

[27! ] . Similar methods can be used to calculate nearest neighbour entanglement at finite temperature [28j |. but the 
conclusion would again be that the threshold is J > T for entanglement. 

We now connect this result to our previous (continuum) investigations |29[ . We can think of the thermal de Broglie 
wavelength for spins as XdB = o-v/ J/T. The condition for entanglement should be that the wavelength is larger than 
the lattice spacing a which leads us to the condition that T < J. Since J = h 2 /2ma 2 , this implies T < h 2 /2ma 2 , which 
is clearly consistent with our previous conclusions (30j . This inequality will now be applied to recent experiments. 



EXPERIMENTAL CONSIDERATIONS 



How would we confirm the existance of spatial entanglement experimentally? To be able to do so we need to 
extract the value of the nearest neighbour two point correlation function (o' 1 + ct^ +1 ). One way of estimating this is 
by interfering independent fluctuating condensates [3l| . Let us now look at the recent remarkable experiment of 
Hadzibabic et al [321 ] and investigate if their experiments reveal any entanglement in two dimensional Bose gases by 
interfering independent condensates. In their experiments the value of the chemical potential is fi/h = lOKHz, while 
the (average) separation between the atoms will be taken to be less than or equal to the healing length a — 0.2/ito 
(as is true for a dilute Bose gas analysed here, see the Appendix). The temperature reached during the experiment 
were 100 — 200nK. They used Rubidium atoms whose mass is 87 atomic units. Putting these numbers together we 
obtain that /i 2 -I- (kT) 2 < J 2 , which as we showed was our criterion for the existence of entanglement. Therefore, 
Bose condensates exhibit field theoretic (continuous variable) entanglement (previously we have only had evidence 
of spin based entanglement in many body systems [23L \3^). Furthermore, we see that the onset of entanglement is 
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very closely related to the BKT transition. The BKT transition was measured to occur at temperatures for which 
the thermal de Broglie wavelength is At = 0.3/im. The healing length, which we can think of as the effective lattice 
spacing, is a = 0.2/im. Entanglement occurs at a < A which agrees with the BKT transition. This can also be obtained 
form the aforementioned field theoretic criterion [2§| that kT < J = h 2 /(2ma), where a is the healing length. What 
is more, we can even estimate the amount of entanglement, as quantified by the single spin reduced entropy. Since 
e± = 1 ± (1 — 2/7rarccos(/i/J)), the (single spin) entanglement is E ss 0.33. The fact that entanglement seems to 
occur at the point where the de Broglie wavelength becomes larger than the healing length, which itself physically 
represents the core size of a vortex, would suggest that entanglement may be linked to the vortex-anti vortex pairing 
that characterises the BKT superfluid phase. This point requires a further more in-depth study. 

CONCLUSIONS 

We have seen that the traditional indicators of order in matter are, superficially speaking, not necessarily strictly 
related to the existence of entanglement between the underlying constituents. However, when the whole situation 
is properly analysed using the information given by a system's Hamiltonian we see that some rigorous relationships 
naturally emerge. We have found conditions for which superfluidity implies entanglement and have then go onto 
demonstrating that some recent experiments have already been sufficiently detailed to witness entanglement in Bosc 
gases. Furthermore, based on very general assumptions, it can be shown that no long range order is possible in one and 
two dimensions, while it is clearly possible to have entanglement under the same circumstances. It is therefore tempting 
to speculate that entanglement is related to low dimensional quantum criticality. The present paper substantiates 
this view with an entanglement analysis of recent cold atom experiments investigating the BKT trantions. Using this 
entanglement for practical purposes would be the next desirable item to investigate in greater detail. 

APPENDIX 
The continuum limit 

We start with the kinetic second quantised Hamiltonian: 

h 2 f 

H = — / &(x)V 2 V(x)dx (9) 
2m J 

where ^(x) represents the annihilation operator for a boson. Suppose that bosons are localised around some lattice 
points (e.g. positions of fixed atoms). Then *&(x) w <fi(x — and ^(x) « cj)*(x — Xi)b\, where 4>{x) are the 
so called Wannier functions. Substituting this into the Hamiltonian leads to 

H = -jJ24bj + Uj2 n i( n i~ 1 ) ( 10 ) 

ij i 

where J = Ti 2 /2m J 4>*{x— Xi)V 2 (j>(x— Xj)dx, U — ATrh 2 s/m J dx\(f>(x)\ 4 and s is the scattering length. If the <fi function 
are well localised than only the nearest neighbour coupling will the relevant, which is what we assumed throughout 
the paper. The on site term is just a constant referring to the total number of atoms. If the density of atoms is low, 
then the nf term can be ignored and we obtain the XX Hamiltonian analysed in the paper. The average energy per 
particle is, at low densities, given by E ss (h 2 /2m)(ps). From here, the healing length, a, is defined by a = l/^fps 
so that the energy is E sa (h 2 /2ma 2 ). This leads to the following witness of entanglement that was mentioned in the 
text: if (H) < (h 2 /2ma 2 ), the system is then entangled 0, [3(3]. This leads to the critical temperature since (H) is 
a function of temperature. Approximately, (H) w ft, 2 /2mA^, where At is the thermal wavelength, so the criterion for 
entanglement is that At > a, which requires the thermal wavelength to be larger than the healing length. 

Two dimensional XX lattice 

True signatures of superfluidity, such as quantized vortices, can really only exist in two dimensions and higher. 
We will see, however, that there will no fundamental changes to our earlier one dimensional conclusions. The two 
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dimensional XX Hamiltonian is given by: 



+ JiK i 4 +1 +< i < j+1 ) (ii) 

We do not know how to diagonalise this Hamiltonian exactly. However, we can do reasonably well by first applying 
generalised (two dimensional) Jordan- Wigner transformation, and then a mean field approximation on the phases. 



The Jordan- Wigner transformation in two dimensions is defined as follows [3J|: a = e %ai >*djj and <r + = e lai ' j d. 



Applying the mean-field on the phases a we obtain the transformed Hamiltonian of the form [341 ]: 

i J 

+ JM,Aj+i -diAj+i) (12) 

which, upon diagonalisation by the two dimensional Fourier Transform, gives 35]: H = A(k)(r}lrik — 5), where the 

eigenvalues are A(k) = \J J\ cos 2 k y + J 2 sin 2 k x . This now presents the (approximately) diagonalised two dimensional 
XX model. Let us now use the energy based entanglement witness [24| to analyse this Hamiltonian. The energy can 
easily be calculated from the partition function and is given by U = — | j A tanh @£ . Applying our energy based 
entanglement witness now reads: \U P \ > (J + J±)/2 implies that the state p is entangled. Let us take the low and 
high temperature limit of the energy formula to draw some conclusions on entanglement. At low T, tanh — ► 1 and 
the (absolute value of) energy is demonstrably larger than the separable bound. The ground state in two dimensions 
is therefore provably entangled. At high T , tanhrr — > x and so \U\ — ► (3(J 2 + Jj_) 2 /8, which implies entanglement for 

temperatures such that kT < g j+j^ i which is very similar to the one dimensional result. 
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